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RESPONSE CHARACTERISTICS OF METEOROLOGICAL ROCKET 


WIND REDUCTION TECHNIQUES 


by 

Alvin J. Miller 

ESSA, Weather Bureau, National Meteorological Center 
Hillcrest Heights, Md, 


INTRODUCTION 

In general, the "mechanical" sources of error inherent in wind 
determinations for meteorological rocket systems can be divided into 3 
categories: (1) sensor response to the wind, (2) radar accuracy of position 
and (3) finite -difference approximation of wind from the position data. 

While Eddy et al. (ref. 2), Malet (ref. 5) and Hyson (ref. 3) have investi- 
gated the first two aspects, the third has received very little attention to 
date. 


Unfortunately, there exists no single procedure for wind reduction, 
mainly because the radar systems and capabilities at each site are so varied. 
Certain stations, for example, are equipped with extremely sophisticated 
radar systems with appropriate Computer equipment so that winds with a 
high degree of resolution result. Other stations must rely on manual re- 
duction of winds from position points depicted on a plotting board at finite 
time intervals (U.S. NASA, ref. 11; Mitchell, ref. 8). 

In the case of tne well equipped stations, the time resolution is, in 
general, such that vertical scales of motion on the order of 1 km are ade- 
quately depicted. Although a serious question exists as to how to smooth 
this detailed wind data or, indeed, if it should be smoothed at all before 
publication in the data books (MRN, ref. 10), we consider this data to be of 
sufficient accuracy for most planned usage (Miller et al., ref. 7). Therefore, 
we will not concern ourselves in this study with these data. In the case of 
the other stations, however, the finite -difference intervals employed in the 
calculations are such that perturbations with "wavelengths" of several 
kilometers are effectively filtered from the data. In view of the results of 
Lettau (ref. 4), Mahoney and Boer (ref. 6), and Cole and Kantor (ref. 1) 
that significant power exists in the wavelength range from 3-4 km, the exact 



degree of smoothing accomplished by the reduction techniques must be known 
if we are to interpret the data correctly. The purpose of this study, then, is 
to investigate the response of the manual wind reduction procedures with 
particular emphasis on the response to small- and medium-scale waves. 

Basically, there are two manual reduction techniques in use at this time, 
height-based and time-based. In the former, (U.S. NASA, ref. 11)* data 
points on the plotting board are interpolated in the vertical so that a position 
and time is depicted at every whole kilometer. Winds are then determined by 
dividing the spatial separation of points 2 kilometers apart by their time 
difference. The wind value is then ascribed to the mid-kilometer level. In 
the time-based technique (Mitchell, ref. 8) position points are depicted at 
30-second intervals for the first 5 minutes from launch, at 60~second intervals 
for the next 15 minutes and at 1Z0 -second intervals thereafter. Winds are 
determined by dividing the position difference by the appropriate time differ- 
ence and the value ascribed to the mean height of the top and bottom points. 
While the time-based method is relatively simple compared to the height- 
based system in that no interpolations are required, it possesses the disad- 
vantage that the vertical distance between points is not constant, since the 
sensor fall rate decreases with decreasing height. 

In this paper analytical expressions are determined for the two methods 
under the assumption that a sinusoidal variation of wind with a fixed vertical 
wavelength exists in the atmosphere. The results are compared and the 
utility of the methods discussed. This analysis is based on a particular fall 
rate profile of the sensor and while this differs somewhat from station to 
station and from instrument to instrument our results should be indicative of 
the problems that arise. 

ANALYTICAL MODEL 

Since our purpose is to evaluate only the wind reduction techniques, we 
assume at the outset that the sensor is perfectly responsive to the winds and 
that the radar positions are exact. 

The expression for the time rate of change of the sensor position (x 
direction only; the results for the y direction would be similar) may then be 
written as: 


ssr <lg , as 

Jhr 

*It has recently come to the attention of the author that the height-based 
technique is to be the sole method employed at all U.S. sites. 
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where x is position, U is the amplitude of the sinusoidal wind variation and 
k is the wave number (k=3ZT. X is the vertical wavelength). 


Clearly, if an analytica^expression for dz/dt as a function of height 
could be determined, then would be analytically defined as a function of 
height. Integration would then yield both x as a function of height and t as a 
function of height. 


For this study, fall rates of Arcasonde-IA instruments with 15 ft. 
diameter parachutes at Wallops Island, Vir ginia (Experimental Inter - 
American Meteorological Rocket Network (EXAMETNET ), ref. 9) were 
plotted and an expression for dz/dt was determined by eye which seemed to 
represent the fall rate adequately. The expression for dz/dt is indicated in 
equation (2) and shown in Figure 1 along with the measured fall rates. 


= (2 ) 
4-c 

where z is in km, A = 2x10”^ km^-sec"^» b = 0. 103 km” ^ , 

While other more complex expressions for dz/dt could be derived that 
would give a somewhat better fit, especially at levels below 30 km, the 
stringent requirement for this study was that finite integrals (non-series 
solutions) exist for both (1) and (2). 

Substitution of (2) into (1) results in the following expressions for x and 
t: 

•*( 2 ) +(2<^(k^))mkK^ 
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Using these exact expressions, the computation of Axt/At by the height- 
based system is relatively straightforward* The values of x and t were 
determined at every whole kilometer and the finite differences calculated as 
described above* However, because (3) and (4) are transcendental, the 
solutions for x and z as functions of t are not readily determined so that for 
the time-based computations a somewhat different procedure had to be 
followed. The values of x and t were determined at 100 -meter intervals 
from 55 km to 25 km and the values linearly interpolated to the necessary 
points in time. While this introduces an uncertainty into the time-based 
calculations it is considered to be small* In all computations the sinusoidal 
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Figure 1. Measured (filled circles) and approximated (solid line) fall rate of Arcasonde 1A instrument 

with 15* diameter parachute at Wallops Island, Va. 



variations were assumed to have unit amplitude 


RESULTS 

Height -based Method 

Illustrated in Fig. 2 are the results of the calculations for the height- 
based method when true wavelengths of 3.5 km and 4.0 km exist in the wind 
profile. The response is about 54% for the former and about 64% for the 
latter; both are essentially independent of height. Of interest is the chopping 
effect that occurs in the 3.5 km wave whenever the true wind maximum occurs 
at a level other than a whole kilometer. This, of course, raises the point 
that for the purposes of our calculations the phase of the sinusoid was held 
fixed. For waves with somewhat different phases, the calculated winds 
would be changed accordingly. 

Figure 3 presents the maximum response for the wind computations, 
when the wind is measured at a level where the true wind is a maximum. 

The response, then, does not include the reduction effect for winds whose 
maxima are at levels other than whole kilometers. For this reason, these 
values should be considered indicative only and are not meant to be used as 
a correction scheme in the wind determinations. 

As might be expected, the response is zero at a wavelength of 2 km and 
increases with increasing wavelength. The effect, then, on vertical scales 
of motion of 3-4 km is to reduce the amplitude by about 36-60%. This, of 
course, has serious repercussions when spectral analyses are attempted on 
this data or when vertical shears are calculated. 

The response at wavelength less than 2 km is somewhat misleading in 
that by this method of wind computation these wavelengths are aliased into 
perturbations of larger scales. Figure 4 presents the results for wave- 
lengths of 1.2 and 1.5 km, selected for illustrative purposes. The aliasing 
of 1.2 km wavelengths into 6 -km waves is clearly evident. Unfortunately, 
information on the reality of these short waves is very scanty at present, 
but indications are (Miller et al. , ref. 7) that they exist with amplitudes of 
the order of 5 m sec”*. Considering however, that all responses for these 
wavelengths are less than 20% the effect of this aliasing does not, in general, 
appear to be significant. 

As a test on the effect of the phase of the wave on the aliasing, compu- 
tations were made for A =1.5 km with the peak of the sinusoid stepped down 
by 100 meter intervals till the next whole kilometer level was reached. No 
significant changes resulted. 
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Figure 3. Maximum response of height-based wind computations as a function of vertical wavelength (AK 







Time -bas ed Method 


/ 

As indicated above, the time -based method involves the use of constant 
time increments that are increased as the fall rate of the sensor decreases. 
While the effect of the selected increments is to allow approximately 2 km 
between data points, it is clear that at the bottom of each time step the 
vertical distance between data points will be less than that at the top. The 
result is a variable smoothing effect in the vertical and we would therefore 
expect the least resolution at the top of each time step and the greatest 
resolution at the bottom. This is borne out by Figure 5 which illustrates 
the computed winds for wavelengths of 1.5 km and 3.5 km. 

In both cases it was assumed that apogee and ejection of the sensor 
occurred within 2 minutes from launch and that 3 minutes of 30 second 
interval data resulted. With the assumption that t=0 at z=55.0 km, the 
heights for the first winds of the second and third time steps are 42.5 km 
and 25.38 km respectively. In the case of the 1.5 km wa,ve, the effect of the 
variable vertical smoothing is clearly seen by the increased resolution at 
the lowest levels, with no pattern evident above about 34 km. With the 
wavelength set equal to 3.5 km, on the other hand, the main disparities 
seem restricted to levels above about 51 km. The effect of the decreased 
resolution at the top of the second step is, however, indicated by the 
relatively low wind value at 42.5 km. Again, the resolution increases with 
decreasing height. 

If apogee were somewhat higher or rocket burn time somewhat longer 
than that assumed above, then the number of data points at the different 
time intervals would be accordingly reduced. As a test on this effect, winds 
were computed in the same way as above, but only 2 minutes of data at 
30-sec. intervals were allowed. The results, also shown in Figure 5, 
indicate a much reduced response between 43 and 47 km compared to the 
previous case. This is due to the increased smoothing interval at these 
heights associated with *he increase in time between data points. 

Because of these disparities and the general increase of resolution with 
decreasing height, no attempt was made to portray a general response curve, 
as was done for the height based system. Clearly, though, the resolving 
power of this method follows a similar pattern with respect to wavelength 
as does the height-based procedure. 
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Figure 5. Winds computed by time-based technique assuming 3 minutes of 
data at 30 second intervals (solid line) and true wind profile (dashed line) 
as function of height. Winds computed assuming 2 minutes of data at 
30 second intervals are indicated by the solid line that merges into the 
stippled line. X represents vertical wavelength of true wind profile. 
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FINAL REMARKS 


By use of analytical expressions for position and time as functions of 
height we have shown that the height-based technique is more consistent in 
portraying medium-scale perturbations than is the time-based method. 

This is accomplished, however, at the expense of the simplicity of the 
time-based system. It is not our intention to argue whether one system is 
effectively better than the other; that is too dependent on the purpose for 
which the data are to be used. Our objective merely is to illustrate the 
differences that result from calculations by the two methods, especially for 
perturbations on the scale of 3-4 km, so that a better understanding of their 
utility might emerge. 

In general, the response of these wind determinations for perturbations 
of 3-4 km wavelength is about 50-60% of true value. This should be con- 
sidered when any attempt at spectral analysis or shear determination is 
made. Both methods alias wavelengths smaller than 2 km with the conse- 
quence that the scale-lengths can not be defined in these data. However 
since both the response to these aliased, waves and the relative power with- 
in them appears to be small, the addition of power to the larger wavelengths 
appears to be negligible. 
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